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Abstract 

A simple algorithm for decoding both errors and erasures of Reed- 
Solomon codes is described. 

1 Introduction 

In this paper, the Gao algorithm modification is given. In the author's 
opinion, the suggested algorithm is the simplest for algebraic codes with short 
lengths for any implementation. 

2 Definitions and notations 



Let us define the (n, k, d) Reed-Solomon code over GF(g) with length 
n = g — 1, number of information symbols designed distance d = n — k+1, 
where q is prime power. 

The message polynomial of the Reed-Solomon code is 

k-l 

M{x) = y^mjx". 

i=0 
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The component Cj of the codeword C{x) is computed as 

d = M{a'), ie[0,n- 1]. 
The received vector is represented as a polynomial 

n—l 71—1 n—1 

R{x) = ^ TjX* = C{x) + E{x) = ^ Qx* + ^ eix\ 

i=0 i=0 i=0 

where C{x) is the codeword, E{x) is the error vector. 

The error vector E{x) has t errors with a set of error positions 
{zi, Z2, • • • , it}- Let us define that Zi = a'^'^ , Z2 = a^^ , . . . , Zt = a** are error 
locations. 

The error locator polynomial is 

t 

W{x) = l[{x-Z,), 
1=1 

where t is the number of errors, Zi is the error location of the error vector 
E{x). 

The error vector E{x) has / erasures with a set of erasure positions S = 
{ill 32, ■ ■ ■ ,ji}- ^1 = C(^^, X2 = a^^ , . . . , Xi = a^' are erasure locations. 
The erasure locator polynomial is 

I 

A{x) = l[{x-X,), 

i=l 

where / is the number of erasures, Xi is the erasure location of the error 
vector E{x). 

The inequality 2t + I < d is well known fT] . 

We construct an interpolating polynomial T(x) such that 

T{a')=ri, ie[0,n-l], 
where deg T{x) < n, and an interpolating polynomial T{x) such that 

r(aO = r,, ie[0,n-l]\S, 
where deg T{x) < n — I. 
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3 Existing algorithms 



We describe here two versions of the Gao algorithm [2], [3l lU [5] . 
The first version is for decoding errors only. Let P{x) = W{x)M{x). The 
key equation is 

{W{x)T{x) = P{x) mod x"" - 1 
degW{x)<^ (1) 
maximize deg W{x) . 

The asymptotic complexity of this algorithm is 0(n(logr;,)^). 
The second version is for decoding both errors and erasures. The key 
equation is 



(W{x)T{x) = P{x) mod 



x"-l 



A{x) 

(degW{x) (2) 
maximize deg . 

The direct computation by this algorithm has complexity O(n^). 
Next, we consider the key equation derivation for the Truong algorithm 
[B] for decoding both errors and erasures. Let 

Q{x) = P{x)A{x) = W{x)M{x)A{x). 

From ([T]) we have 

W{x){{T{x)A{x)) = (P(x)A(x)) mod — 1 
and the key equation is 

{W{x) ((T(x)A(x)) = Q{x) mod x" - 1 
degW{x)<^ (3) 
maximize deg W{x). 

The asymptotic complexity of this algorithm coincides with the complex- 
ity of decoding algorithms [21 [3l HJ Ej . 
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4 Suggested algorithm 



We introduce the following lemma. 

Lemma: 




Proof: From Newton's interpolation formula we obtain 



T{x) 



X' 



- 1 



u{x) + r(x) 



A(x) 



where U {x) is a polynomial. 

From (111) and the lemma we get a new key equation 



W{x)T{x) = P{x) mod 
deg W{x) < 
maximize deg W{x) . 



(4) 



The description of the three algorithms for decoding both errors and 
erasures is in table 1. 

5 Conclusion 

The suggested algorithm has replaced the computation using Newton's in- 
terpolation formula by the fast computation of the discrete Fourier transform. 
The algorithm complexity is less than the Truong algorithm [6] complexity 
because the suggested algorithm does not contain some of the intermediate 
computations. 
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Table 1 

Algorithms for decoding both errors and erasures 



Step 


Gao's algorithm 


Truong's algorithm 


Suggested algorithm 



1 

2a 
2b 
3 


x"- 1 
Kix) 

{w{x)T{x) = P{x) 

J A{x) 

I deg W{x) < '^-^-^ 
maximize deg W{x) 

M(x) = —-4 


A{x) 
T{x) 

T{x)A(x) 

(W(x){(T{x)A(x))=Q(x) 
1 moQ X — 1 
1 degW{x) < '^-^-^ 
maximize deg W{x) 

Mix) = ^^""^ 


nx) 
- 1 

Mx) 

(w(x)T(x) = P{x) 

J A(a;) 

I degW{x) < ^-^-^ 
maximize deg W{x) 

Mix) = , ; 

^ ' W{x) 


Complexity 




0{n{\ognf) 


0{n{\ognf) 
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